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Abstract : The generalized logarithmic series distribution (GLSD) adds an extra parameter to the 
usual logarithmic series distribution and was introduced by Jain and Gupta (1973). This distribution 
has found applications in various fields. The estimation of parameters of generalized logarithmic 
series distribution was studied by the methods of maximum likelihood, moments, minimum chi-square 
and weighted discrepancies. The GLSD was fitted to different leaf spot grades in four varieties of 
mulberry namely Ichinose, Gosherami, Rokokuvoso and Kokuso-20 and it was observed that all the 
estimation techniques perform well in estimating the parameters of generalized logarithmic series 
distribution but with varying degree of non-significance. 
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Introduction : 

The generalized logarithmic series distribution 
(GLSD) characterized by two parameters a and |3 was 
first obtained by Jain and Gupta (1973) and its probability 
function is given by: 


(1) 

o< <i, | |<i 

Since GLSD is generalization of logarithmic series 
distribution (LSD). GLSD will reduce to logarithmic series 
distribution (LSD) when taking P= 1 . GLSD is a member 
of Consul and Shenton (1972) family of Lagrangian 
probability distribution and also Gupta (1974) modified 


power series distribution (MPSD). The same distribution 
has been obtained by many more authors. The distribution 
has found applications in various fields. Jain and Gupta 
(1973) applied it to the William’s data on number of papers 
by entomologists, Rao (1981) applied it to the study of 
correlation between two types of children in a family, 
and Hanseen and Willenkens (1990) used it in the risk 
theory in a problem related to the total claim size upto 
time t. The estimation of GLSD has been studied by many 
researchers, whereas Gupta (1974) and Jani (1977) 
examined its minimum variance unbiased (MVU) 
estimation. Mishra (1979) discussed its maximum 
likelihood (ML) estimation, Jani and Shah (1979) 
discussed the method of moments of its estimation. 

Since GLSD is generalization of logarithmic series 
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distribution (LSD). GLSD will reduce to logarithmic series 
distribution (LSD) when taking P= 1 . 

Put P=1 


p(X = x) = 



£ 

IxiVT 


*(!- )° 


i (x-i): x 
log(l- )x(x-l) 


LogL = Nlog + N x log + S V f x log ( x-j)+( -l)Nxlog(l- LZfilogX; 

X=1 j=l 

The two likelihood functions can be obtained as 
aLogL N Nx ( - 1) Nx ft 

a = ~(i- ) + “ (i- ) =0 

aLogL _ - Nx | k x-l xf x 0 

a x=i j=i x-j 


( 2 ) 

(3) 


So one parameter LSD is a particular case of GLSD 

1 X 

xlog(l - ) 

Materials and Methods : 

For the present study mulberry data set was utilized. 
Four different varieties of mulberry namely Ichinose, 
Goshoerami, Rokokuyoso and Kokuso having different 
leaf spot disease intensity were chosen for the study. 
Three trees of each variety were selected at random, 
from each tree, three branches were selected at random 
and then from each branch the spots of the disease were 
recorded from all the leaves. The leaves with no spot 
were referred to as disease free and named as grade 
zero (0 grade). The leaves having 1 to 5, 6 to 10, 1 1 to 15, 
16 to 20 and more than 20 spots were total graded as 1, 
2, 3, 4 and 5 grade, respectively. The number of leaves in 
three branches per tree falling in different grades was 
recorded. Similarly for all the four varieties, the number 
of leaves in three branches falling under different grades 
were recorded. All the leaves of the three randomly 
selected branches per plant were observed for different 
leaf spot grading. 

Estimation of parameters : 

Maximum likelihood estimates : 

Consider a random sample of size N taken from the 
GLSD and let the observed frequencies be 
k 

f x ; x = 1,2. ..k so that £ f x = N where k is the largest 

X=1 

of the observed values having non-zero frequencies. The 
likelihood equation of the GLSD can be written as 

N Ns nn(x-j)V ) ( ~ 1)Nx 

L = 

rpx,!/* 

i=l 

The log likelihood function is given as 


where x is the sample mean. From equation (3), 
we get 


i _ 

X 

Putting this equ. in (3), we get 



k x-l 

Z Z 


x=lj=l 



= 0 


(4) 

(5) 


Which can be solved for a and p the M.L. 
estimators of a and p, respectively 


Estimates from moments : 

The mean, the variance and the recurrence relation 
for higher moments is given by 


1= ^ 



r+1 


iLi 



3 r-1 


(6) 

(7) 

( 8 ) 


Since p 2 < 0 for p0 < 1 - (l - 0)/ a0 , we will 
obtain the moment estimators of the parameters of the 
GLSD for the sample space 1 - (l - 0)/ a0 < P0 <1 
From eq. (6) and (7), we have 
2 2 - Q (l - )=o (9) 

where, 

Q= 3 /( 2 +2 ) (10) 

The equ. (9) can be solved for 0 by using the method 
of iterations. To obtain the initial value of 0 expanding 

a 2 0 2 -Q(l-0) into power series expansion and 
neglecting 0 3 and terms higher than that we get (9) as 
2 - 12(Q -lXl - ) = 0 
which gives, for 0>O 

= 2(VQ'(Q' + 1))-Q') (11) 


where, Q' = 3(Q - l) and Q is given by equ. (10) 
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To get 9*, the moment estimator of 9, p and p, are 
to replaced by their respective estimates sample mean 

x and sample variance S 2 of the observed data. 

Using equ. (6) and replacing p by x.we get 

’= 1 / *- 7x (12) 

Minimum chi-square for GLSD: 

Since 


k. Obviously, k is the largest of the observations. 

Let n = if x . 

x=l 

The corresponding relative frequencies are given 
by 

n x =^-; X = 9,1,2,... k (17) 

The log likelihood function can be written as 

logL= ZNn x logp x ( , ) (18) 

X=1 


2 _ y k -Px) 2 
x=l Px 

is approximately distributed as chi-square. 
Differentiating it with respect to a and p, we obtain 


I k -Pxli+— 1^-iogPx =<> 

Px J 8 


Z K -Pxfl + — T“ lo §Px =0 
x=l l Px J 8 


and (14) we get 


k / N n 

I («x-Px) 1 + — 
X=1 V Px 


Nx ( -l)Nx 


(l- jlogll 


k 

z 

X— 1 


Z(n x -P x ) 1 + — 

P 


— , , k x— l xn 

Nxlog(l- )+ I Z ■ 


x=lj=l( X-j 

The resulting equ. (15) and (16) are known as 
minimum chi-square equations. 

Weighted discrepancies (WD) method : 

Let f denote the observed frequencies x=9, 1, 2, 


The likelihood equations are 


k 8 

Zn x — logp x =0 

x=l O 

k d 

I n x — logp x =0 

x=l o 


(19) 

(29) 


(13) 

k 

Again as Lp x = i 

k a 

(21) 

(14) 

Z Px— logPx =0 

X— 1 O 

(22) 

to (13) 

k a 

Z Px— !°gPx =0 

x=l o 

(23) 


Subtracting (22) from (19) and (23) from (26), we 
get 

(15) 

k , x d 

(24) 

IK-Pxj— ^ logPx =0 

x=l O 

(16) 

k / \ Q 

Z( n x -p x )— Iogp x - 0 

Y-1 O 



(25) 


Substituting the corresponding expressions of the 
derivatives to (24) and (25), we get 


S( n x-Px 


N 

Nx ( -l)Nx 

(1- >og(l- 

r (i- ) J 


(26) 


Table 1 : Comparison of observed frequencies with expected frequencies of GLSD for leaf spot grades in Ichinose variety of mulberry 

Expected frequency of GLSD 

Leaf spot grades Observed frequency Methods of estimation 


ML Moments MC WD 


0 

61 

64.70 

62.80 

64.00 

62.50 

1 

18 

16.10 

15.30 

13.15 

15.50 

2 

1 1 

9.10 

10.55 

10.45 

10.50 

3 

3 

2.00 

3.30 

3.50 

3.40 

4 

1 

1.10 

1.10 

1.70 

1.20 

5 

0 

1.00 

0.95 

1.20 

0.90 

Total 

94 

94.00 

94.00 

94.00 

94.00 

Parameters (Estimate) 






a 


0.871 

0.867 

0.862 

0.853 

P 


0.961 

0.921 

0.932 

0.915 

l 2 


1.645 

1.305 

2.402 

1.145 

p-value 


0.896 

0.934 

0.791 

0.950 
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iK-p^niMi- (27) 

which has referred to Kemp (1986) as an equation 
from minimum discrimination information and ML 
estimation and called as weighted discrepancies 
estimation method. 

Application : 

The GLSD has been found useful in research areas 
where the logarithmic series distribution has been applied. 
Jani and Shah (1979) used the GLSD to model the count 
of European red mites on apple leaves. The different 
methods of estimation for estimating the parameters of 
generalized logarithmic series distribution were compared 
by chi-square goodness of fit test. 

In the Table 1 the estimates of parameters a and |3 
in different methods of estimation are ML (0.87 1 , 0.96 1), 
Moments (0.867, 0.921), MC (0.862, 0.932) and WD 
(0.853, 0.915). The p-value of all the methods are non- 
significant and in agreement with the observed 
frequencies. We observed that all the estimation 
techniques, the ML, the moment, the MC and the WD 
method perform well in estimating the GLSD parameters. 

In the Table 2 the estimates of parameters d and |3 
in different methods of estimation are ML (0.861, 0.952), 
Moments (0.845, 0.931), MC (0.855, 0.940) and WD 
(0.832, 0.925). The p-value of all the methods are non- 
significant and in agreement with the observed 
frequencies. We observed that all the estimation 
techniques, the ML, the moment, the MC and the WD 
method perform well in estimating the GLSD parameters. 

In the Table 3 the estimates of parameters d and (1 
in different methods of estimation are ML (0.870, 0.955), 
Moments (0.810, 0.915), MC (0.865, 0.950) and WD 
(0.862, 0.941). The p-value of all the methods are non- 
significant and in agreement with the observed 
frequencies. We observed that all the estimation 
techniques, the ML, the moment, the MC and the WD 
method perform well in estimating the GLSD parameters. 

In the Table 4 the estimates of parameters a and (3 
in different methods of estimation are ML (0.775, 0.915), 
Moments (0.785, 0.852), MC (0.720, 0.920) and WD 
(0.753, 0.849). The p-value of all the methods are non- 


significant and in agreement with the observed 
frequencies. We observed that all the estimation 
techniques, the ML, the moment, the MC and the WD 
method perform well in estimating the GLSD parameters. 
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